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Abstract 

Q^, In this article, we study a second-order expansion for the effect induced on a large 

quantum particle which undergoes a single scattering with a low-mass particle via a re- 
pulsive point interaction. We give an approximation with third-order error in A to the 



map G Tr2[(/ ® p)S\{G ® I)Sx], where G € B{L'^{W^)) is a heavy-particle observable, 



m 



p G i?i(M") is the density matrix corresponding to the state of the light particle, A — ^ 
is the mass ratio of the light particle to the heavy particle, Sx € i?(L^(M") (S) L^(M")) 
, is the scattering matrix between the two particles due to a repulsive point interaction, 

\^ [ and the trace is over the light-particle Hilbert space. The third-order error is bounded in 

operator norm for dimensions one and three using a weighted operator norm on G. 

O ■ 1 Introduction 

oo ■ 

In theoretical physics, many derivations of decoherence models begin with an analysis of the 
effect on a test particle of a scattering with a single particle from a background gas P [6l [8]. 
^ ■ A regime that the theorists have studied and which has generated interest in experimental 
physics [7] is when the test particle is much more massive than a single particle from the 
gas. Mathematical progress towards justifying the scattering assumption made in the physical 
literature in the regime where a test particle interacts with particles of comparatively low mass 
can be found in |ll[3],|5]. In this article, we study a scattering map expressing the effect induced 
on a test particle of mass M by an interaction with a particle of mass m = AM, A -C 1. The 
force interaction between the test particle and the gas particle is taken as a repulsive point 
potential. 

We work towards bounding the error e{G, A) in operator norm for G G -B(L^(R")), tt, = 1, 3 
of a second order approximation: 

Tr2[(/ ® p)SliG ® /)Sa] =G + XMi{G) + X^N^G) + e{G, A), (1.1) 

where p G Bi{L'^{R'')) is a density matrix (i.e. p > and Tr[p] = 1), G G B{L'^{W)), 
Sx G i?(L^(R"') ® L^(M"')) is the unitary scattering operator for a point interaction, and the 
partial trace is over the second component of the Hilbert space L^(M") ® L^(R"'). Mi and M; 
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are linear maps acting on a dense subspace of i?(L^(R")) {M2 is unbounded). Our main result 
is that there exists a c > such that for all p, G, and < A 

||e(G,A)|| <cA^|| licil 

where || ■ is a weighted operator norm of the form 

\\GLn = IIGII + |||X|G|| + 
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0<i,j<d ei+e2<3 

and II ■ Wwtn is a weighted trace norm which will depend on the dimension. In the above, X and 
P are the vector of position and momentum operators respectively: {Xjf){x) = Xjf{x) and 
{Pjf){x) = i{-^f){x). Expressions of the type A*GB for unbounded operators A and B are 
identified with the kernel of the densely defined quadric form F{i(ji; 'i/j2) = {AiIji\GBiIj2) in the 
case that F is bounded. 

The scattering operator is defined as Sa = where 

= s-lime**-^*°'e-**-^'='" (1.2) 

t— >±oo 

are the Moller wave operators, and Hkm is the kinetic Hamiltonian and is the standard self- 
adjoint extension of the sum of the Laplacians —^A^eavy — ^Ang^t, while the total Hamilto- 
nian Htot includes an additional repulsive point interaction between the particles. The definition 
of Htot is a little tricky for tj, > 1 since, in analogy to the Hamiltonian for a particle in a point 
potential [2j, it can not be defined as a perturbation of Hkm even in the sense of a quadratic 
form. Rather, it is defined as a self-adjoint extension of —^Aheavy — ^Aught with a special 
boundary condition. Going to center of mass coordinates, we can write 

1 1 A _ 1 A M + m 

— Aneavy + l^^l^9ht - + M) ^™ 

SO that the special boundary condition will be placed on the displacement coordinate corre- 
sponding to Adis and follows in analogy with that a single particle in a point potential as 
discussed in [2]. This also allows us to write down expressions for Sa. Non-trivial point po- 
tentials in dimensions > 3 do not exist and the main result of our analysis is restricted to 
dimensions one and three. 

The first and second order expressions Mi{G) and M2{G) respectively have the form 

Mi(G') = i[V,, G] and M2{G) = 1^2 + \{A, P}, G\ + ^{G) - ^ip{I)G - ^Gip{I), (1.3) 

where Vi, V2, and {A)j for j = 1, . . . ,n are bounded real-valued functions of the operator X, 
and (/9 is a completely positive map admitting a Kraus decomposition: 

^{G) = Y, I dkml^Gm^^^, (1.4) 

with the nij ^'s being bounded multiplication operators in the X-basis. Notice that terms 
in ( ll.3p are reminiscent of the form of a Lindblad generator [10]. In [1] the results of this 
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article are applied to the convergence of a quantum dynamical semigroup to a limiting form 
with generator including the terms (11.31) . 

The explicit forms for Vi, V2 , A, and f are: 

Vi = CnS-^ [ dk\k\-^ [ dvidv2p{vi,V2)e'^'^'^'-'^'\ (1.5) 

V2 = CnS;;^ f dk\k\-^ [ dv^dV2iVi + V2)VTp{VuV2)e''^^'^'-'^'\ (1.6) 

A= CnS-' [ dk\k\-^ [ dvidv2VTp{vi,V2)e'^^'^'-'^^\ (1.7) 

^{G) = els-' [ dk\k\-' [ yv,dv2p{vi,V2)e''^'^-'''+'^^Ge-''^^-^'+'^\ (1.8) 

where is surface area of a unit ball in M", c„ is a constant arising form the scattering operator 
Sa, p{ki, ^2) is the integral kernel of p, and Vt is the gradient of weak derivatives in the diagonal 
direction which is formally (Vtp(^i, ^2))^ = lim/i^o (p(^i + hej,k2 + hej) — p{ki,k2)). 
The integral kernel p{ki,k2) is well defined since p is traceclass and hence Hilbert- Schmidt. 
In dimension one, the integrals J^^^|^|^2|=fc replaced by discrete sums. In dimension two, 
V2 has an additional term due to the logarithm in (14.41) which we did not write down in 
the expression for V2 above. The multiplication operators nij ^ are defined as rrij ^(X) = 
CnS~'^ J\^l=\k\ fi^^) (,-iX(-v+k)^ where p = Y^j f^j\fj){fj\ is the diagonalized form of p. Vi, 
V2, A, and ip are bounded under certain norm restrictions on p, since, for example, ||Vi|| < 

C„|||P|"-2p||, and IIV^II = C„|||P|"-2p|p|n-2||^^ 

With the center of mass coordinate at the origin, the scattering operator S (neglecting the 
index A until it is explained) acts identically on the center-of-mass component of L^(]R^") = 
L2(R") ® L2(R") as 

S = I + hm^{S{k)0\(l>){ct>\), (1.9) 

where the right copy of L^(M") corresponds to the displacement variable and is decomposed in 
the momentum basis into a radial and an angular component as L''{M.~^ ,r"'-^dr) L''{dBi{0)), 
S{k) acts as a multiplication operator on the L^(M+) component, and = (snJ^'^ldBiio), is 
the normalized indicator function over the whole surface dBi{0). We call ^(A;) the scattering 
coefficient, and it has the form 

Dim — 1 : Dim — 2 : Dim — 3 : 

SM = T^^ Si{k) = - 5,(A:) = — (1.10) 

k + i^a + 7 + ln(f) + «f ' ' l-^ + ik 

where a is a resonance parameter defined for the one-dimensional case, / is the scattering length 
in the two- and three-dimensional cases and 7 ~ .57721 is the Euler-Mascheroni constant. In 
the one-dimensional case a scattering length / is sometimes defined as the negative inverse of 
the resonance parameter a = j^, where c is the coupling constant of the interaction and p is 
the relative mass ^^ttt = ^ttt- However, this contrasts with the two- and three-dimensional 

m+M 1+A ' 

cases where the scattering length is proportional to the strength of the interaction. In the 
context of this article, where the point interaction is between a light and an heavy particle. 
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we parameterize the resonance parameter as a = j^c^q in the one- dimensional case and the 
scattering length as / = j^ki in the two- and three-dimensional cases for some fixed ao and 
/. This corresponds to holding the strength of the interaction fixed. Thus Sa and Sx will be 
indexed by A for the remainder of the article. 

There are two main obstacles in attempting to find a bound for the error e{G, A) from (11.11) . 
The first obstacle is to find helpful expressions to facilitate making a Taylor expansion in A of 
Tr2[(/ (S) p)Sl{G ® /)Sa]. Writing = Sa - /, then 

Tr2[(/ ® p)Sl{G ® /)Sa] = G + Tr2[(/ ® p)Al]G + GTi^Hl ® p) Aa] + Tr2[(/ ® p) A* (G ® /) Aa] , 

and it turns out to be natural at all points of the analysis to approach the terms on the right 
individually. Propositions 12.21 and 12.31 are directed towards finding expressions for 

Tr2[(/ ® p)A*] and Tr2[(/ ® p)A*(G' ® /)Aa] (1.11) 

respectively (since Tr2[(/ ® p)Aa] is merely the adjoint of Tr2[(/ p) A^]). The expressions we 
find in Propositions 12.21 and 12.31 are of the form 

Tr2[(/®p)A*]G= / dk [ daUIJi^G, and (1.12) 
Jm." Jso„ 

Tt,[{I®p)AI{G®I)A,]= [ dk [ da,da,U*^ hi Gh^,^^^^,U^,^^^^„ (1.13) 

for some unitaries U*. .Ur x and some bounded operators or ^, h*. ^ which are functions 
of the vector of momentum operators P. In general, we will have the problem that 

dk da \\g^^^ J = oo, and dk daida2\\h^^^^J\\h^^^^J = oc, 

so the integrals of operators only have strong convergence. Propositions 13. II and 13.31 make sense 
of the integrals of operators such as fll.l2p and (I1.13P that arise and give operator norm bounds 
for the limits. The basic pattern in the proof of Propositions 13.11 and 13.31 is an application of 
the simple inequalities in Propositions IA.2I and IA.3I in addition to intertwining relations that 
we have between the multiplication operators and the unitaries appearing in fll.l2p and (11.130 . 

Bounding the third order error of the expansions in A of the strongly convergent inte- 
grals (I1.12P and (I1.13P brings up the second major obstacle. We will need to bound certain 
strongly convergent integrals for all A in a neighborhood of zero. For small A there will be 
unbounded expressions arising from the scattering coefficient Sx{k) that will have contrasting 
properties between the one- and three-dimensional cases. For example, In the limit A — > 0, 
jS\{k) becomes increasingly peaked in absolute value at /c ~ in the one- dimensional case. 
For the three-dimensional case, jS\{k) becomes increasingly peaked at A; = oo. A difficulty 
with the two-dimensional case is the presence of the natural logarithm in the expression for 
Sx{k) and the fact that jS \ j^(k) is not peaked at a fixed point as A varies. The peak point 
does tend towards A; ~ as A — 0, but it is unknown how to attain the necessary inequalities 
in this case. 

This article is organized as follows. Section [2] is concerned with proving Propositions 12.21 
and 12.31 which give expressions for Tr[(/ ® p)A^] and Tr[(J ® p)A^(G' ® /)Aa]. In Section [3] 
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we prove Propositions 13.11 and 13.31 which give the primary tools for bounding the integrals of 
operators which will arise in bounding the error term e{G, A) of our expansion (11.11) . Section H] 
contains the proof of Theorem 14 . 21 which is the main result of the article. This involves expand- 
ing the expressions in Propositions 12.21 and 12 .31 that we found in Section [2] in A and bounding the 
error. The difficult parts of the proof are characterized by using the Propositions 13. II and 13.31 to 
translate unbounded expressions arising from the expansion of the scattering coefficient Sx into 
conditions on G and p through the weighted norms HGHiun and ||p||u,jn being finite. Sections [2] 
and [3] apply to dimensions one through three (all dimensions where non-trivial point potentials 
exist), while Section H] does not treat dimension two. 



2 Finding useful expressions for a single scattering 

In this section, we will find expressions for Tr2[pA^] and Tr2[pA^G'A;v]. For notational conve- 
nience, we will begin identifying p with p and G®I with G. Finding formulas for Tr2[pA^], 
Tr2[pA^GAA] begins with writing Aa = Sa — / in a convenient way. Let f,g & L^(M" x M"), 
where the first and second component of x correspond to the displacement and the center 
of mass coordinate, then 

{g\Axf)= [ dK,mrdk^^{[ dhgCki,Kcm)){ [ dk fik, Kcm)) ■ (2.1) 

The above formula gives a quadratic form representation of Aa that involves integrating over a 
surface of 3n — 1 degrees of freedom rather than An, since it acts identically over the center-of- 
mass component of the Hilbert space and conserves energy for the complementary displacement 
coordinate. The integral kernel for Aa in center-of-mass momentum coordinates can be formally 
expressed as 

A\{kdis,l, Kcm,l'i kdis,2, Kcm,2) = TT T^^'^l | | ~ | ^dis,2 1 )^(-^cm,l ~ -^cm,2) 

I '^dis,l\ 

However, for instance, this does not work directly towards finding even a formal expression for 
{Tr2[pAl]G){K^,K2) = j dk, dh dK p{k,, k^) A^ih, K^h, K) G{K, K^), (2.2) 

where we have written down a formal equation between integral kernel entries Tr2[pA^G], G, 
and A^ using momentum coordinates corresponding to the heavy particle and the light particle. 
In finding an expression for (12.21) . it would be natural to have K2 as a parameterizing variable 
since the expression above is just multiplication of G from the left by Tr2[pA^]. 

For A = the center of mass coordinates are Xcm = + Xd = x — X, 

where x and X are the position vectors of the particle with mass m and M. The corresponding 
momentum coordinates are kd = x+a^ ^ TTa = k + K. The proposition below 

gives two quadratic form representations of Aa using different parameterisations of the inte- 
gration in (12.11) . (12.31) is directed towards finding an expression for Tr2[pA^] and (12.41) is for 
Tr2[pA^GAA]. The proof of the following proposition requires changes of integration. 

Proposition 2.1 (Quadratic form representations of Aa). Let f,g& L^(]R" x M"), then 
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1. First Quadratic Form Representation 



{g\Axf)= / dkdK / daSx{\k\) 

Jr" J SOn 

g{k + X{K + ak), K + {a - I)k) f{ak + X{K + ak), K), (2.3) 
2. Second Quadratic Form Representation 

{g\Ai,f)= I dK^dh I dadet{I + Xa)-'Sxi\j^{h-XK2)\) 

9{ki, K2 + - Ai?2)) fih + ^^ih - XK2), K2), (2.4) 

1 + Aa 1 + Act 

where the total Haar measure on SOn is normalized to bel ( and for dimension one, the integral 
over SOn is replaced by a sum over {+, —}). 

The proofs of Propositions 12 . 21 and 12 . 31 work by using the spectral decomposition of p, special 
cases of G, etc. so that the quadratic form representations (12. 3p and (12.41) of can be applied. 
Defining = e'^'''^ , recall that acts in the momentum basis as a shift: {Tj:f){p) = f{p — k). 

Proposition 2.2. Let p have continuous integral operator elements in momentum representa- 
tion. Tr2 [pAl\ has the integral form 

Bl= I dkj rf^r,^r;^^P,v,A^A(|fc|), (2.5) 

JR" J SOn 

where T3 is a translation by a in the momentum P basis and pj:^^ is a multiplication operator: 

Pk,a,x = P((l + A)fc + AF, (a + A)^+ AF). 
Proof. The following equality holds: 

Tr^bA*] = Tr2E/3,|/,)(/,|A*] = /?,(id ® (/,|)Al(id ® |/,)), 

j j 

where the infinite sum on the right converges absolutely in the operator norm. If we take a 
partial sum = Y.'jli (^j\fj){fj\^ then using ([23D, 

^(«;|(id® (/,|)A*(id® |/,))i.) = J2 dKidK2 / dkSxm / da 

fj{ak + X{Ki + ak)) w{Ki) f,{k + A(i?i + ak)) v^K^ + (a - I)k). 
This has the form where [■] is given by 

[ dkSxik) I daT*^Tf:pU{l + X)k + XK,{a + X)k + XK). 

JR" JsOn 

This converges in operator norm to the expression given by (12. 5p . since pm — ^ p in the trace 
norm and by the bound given in Corollary 13. 2[ 

□ 
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Tr2[pAA] has a similar integral representation by taking the adjoint. Now we will delve 
into the form of Tr2[pA^G'AA]. In the following, the operator Da acts on / G L'^{MP') as 
{DAf){k) = I det{A)\^f{Ak) for a element A G G'L„(M). 

Proposition 2.3. Let Ylj l^j\fj){fj\ spectral decomposition of p. Tr2[pA^G'AA] can be 

written in the form 

f f - k — XP 

B,(G) =J2 dk 5.(1^—^1) 
■ Jm" Jso„xso„ i + A 



/^det(l + Aa)-5/^.(^+ ^(^_ AF)). 



k — \ P 

G^A(|^:py|)m,,>„.f/,V,„ (2.6) 

where Ur , = rtD i+xa T r,. T„h, rr, and D i+\ act on the momentum basis and m.r . is a 
function of the momentum operator P of the form 

- / 

T 

Proof. Equation fl2.4l) tells us how A^ acts as a quadratic form. In order to use (12.41) . we will 
look at {v\TT2[pA\GAx]w) in the special case where G = G^I = \y) {y\ ®/ is a one-dimensional 
projection tensored with the identity over the light-particle Hilbert space. Formally, this allows 
us to write 

{v\TT2[pAlGA),]w) = ^ ® /,-|A^|2/ ® ® (j)i\Ax\w ^ fj), 

j I 

where (0m) is some orthonormal basis over the light-particle Hilbert space allowing a repre- 
sentation of the identity operator as a sum of one- dimensional projections , and the spectral 
decomposition of p has been used. Once ( 12. 4p has been apphed, we build up to an expres- 
sion (12.61) . taking care with respect to the limits involved. By Corollary 13. 41 the expression (12. 6p 
defines a bounded completely positive map (c.p.m.). Since Tr2[pA^GAA] defines a c.p.m. and 
agrees with (12. 6p for one-dimensional orthogonal projections, it follows that the two expressions 
are equal on B{L'^{W^)). This follows because c.p.m.'s are strongly continuous and the span of 
one-dimensional orthogonal projections is strongly dense. 

The following holds, where the right-hand side converges in the operator norm: 

Tr2[pAlGA,] = ^/?,(id® (/,|)A*G'A,(id® I/,)). 

j 

For G = \y){y\, (id {fj\)AlGAx{id (g) \ fj)) = (pyj{T), where Lpyj is the completely positive 
map such that for H G BiH) 

Vy,{H) = (id ® {f,\)Al{\y){y\ ® //)A,(id ® \f,)). 

Since Lpyj is completely positive, '^yji^JLi I0i)(0«l) converges strongly to ipyj{I). (Pyj{I) is 
determined by its expectations {v\ipy{I)v), and moreover 



N 

{v\^yj{I)v) = lim {vl^yjiY] \(j)m){(pm\)v) 

m=l 

^ f - - - 

\im y2{(j)^m\vv,j,y){v^,j,y\^„,) = \\vyj,yf = / dk J ^y{k) J ^y{k) , (2.7) 



N^OO 

m=l 



7 



where v^j^y is defined as the vector v^j^y = {{y\ ^id)Ax{\v) ® Using (12.41) . {(pm\vv,j,y) can 
be expressed as 



V, 



,j^y)= [ dK [ daSx{\ \ (fc- Ai?)|) det(/ + Aa)-^ 
J JsOn / + Act 



4>m{k) y{K+ ^-^ik - \K)) f, {k + ^-^{k - XK) v{K). (2.8) 
1 + A(j 1 + A<y 

By (12. 7p . we can evaluate Tr2[(|/j) (/j|)A^(|?/) (y| ® I)-A\\ = {'v\'^y{I)v) through expression 
dkvyj^y{k)vvj^y{k). Through (12.81) we have an a.e. defined expression for the values Vyj^y{k). 
Now, writing down dkvyj^y{k)vyj^y{k) using the expression for Vyj^y{k), the result can be 
viewed as an integral of operators acting from the left and the right on \y){y\, followed by an 
evaluation {v\{-)v). Using the intertwining relation: 

miP)TpDi+x. T r = TpDi+x. T tmiP- ^~ ^ (k- XP)), 

for a function m{P) of the momentum operators P and the fact that = a ^Jr^^^^ is an 
isometry for < A < 1, the expression can be written: 

r _ r _ ^ _ \p 

{v\ipyjil)v) = {v\ dk dai da2 [Ul ^ m* g Jx[\ | ) 

k — XP 

{\y){y\)Sx{\ I ) "^jAa^ ,A f^fc,.^ ,a] I ^) • 

So fy,j{I) = Tr2[(|/j)(/j|)A^(|t;)(t>|)AA] agrees with the expression (12.61) for a fixed j and for 
G = \v){v\ forallf, and hence by our observation at the beginning of the proof, TY2[{\fj){fj\)A\GAx] 
is equal to the expression (12. 6p for a single fixed j and all G G i?(L^(]R")). However, if we take 
the limit m — oo for = XljLi l^j\fj){fj\y then the expression (12.61) converges in the operator 
norm and Tr2 [pm A^GA^] converges to Tr2[pA^GA;y]. Hence we have equality for all trace class 
P- 

□ 

Through the formula Tr2[/9S^GSA] = G+ B*G+GB + B{G), it is clear that B* + I3 = -B(/) 
by plugging in G = /. However, it is not at all obvious that this equality takes place through 
the expressions (12.51) and (12. 6p for B* and B(J), respectively, since the operators U^:^^ appear 

only in form for B(/). 

It is convenient to notice the intertwining relation h{k — XP) Uj:^^ = U^^^ h{j^^k — XP)) . 
Let g E L^(]R"'), then g = 'B{I)g can be written: 



Kp) = J2 m rn*-^ U^^^^ ,) (p) 

■ JR" JSOnXSOn 



\Sx 



I 



ilj^ik - Xp)\) {Ul^^^m^^,^^^^,^^^^^^^ (2.9) 



where we have intertwined Ui from the left to the right, and 

fc,cri,A 



i^Ux "^k..A Ul^J (P) = det(/ + AaO-^ /, {k + " Ap)) . 
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iUUxU^,,..,x9m = det (7^)^ det C-^)"' 9{p+ (-1 " + Xa,)-\k - Xp)). 

Making the change of variables j^^{k — Xp) —>■ k, the resulting expression has only angular 
dependance of a"20"f ^ = a, and integrating out the other angular degrees of freedom yields B. 

3 Bounding integrals of non- commuting operators 

Now we move on to proving Propositions 13.11 and 13.31 below which are proved in much greater 
generality than needed for this section, but they will serve as the principle tools in Section |H To 
state these propositions we will need to generalize the concept of a multiplication operator. Let 
TCi, TC2 be Hilbert spaces. Given a bounded function M : M" — B(Hi,H2) we can construct an 
element M G i3(L2(M") ® Hi, ^^(M") ® 7^2) using the equivalence L'^(R") ®ni = L'^{W, Hi), 
where for f e L'^{W)®ni 

M(f)(f) = M(f)f(f). 
We will call these multiplication operators. 

Proposition 3.1. Define B : L'^iW) (^Hi^ L2(M") ® s.t. 

B= dk darlr^jqj:,^, (3.1) 

where q^^ is a multiplication operator in the P basis of the form: 

%a = '^k,a(^)v{^l,cTk + y^P, X2,ak + JaP), 

where r]{ki,k2) is continuous and defines a trace class integral operator on L^(]R"'), 

a^r, xi^o-, X2,o-, Ycr G M„(M), andnj:^ G i3(L^(]R"') (g)7ii, L^(M") (g) 7^2) is a multiplication operator. 

Let 

I det(xi^^ + y<^(a^ - det(x2,<^ + y<x(a^ - det(xi^^)|, and \ det(x2,<^)| 

be uniformly bounded from below by - for some c > 0. Finally, let the family of maps n^^{K) G 
-8(7^1,7^2) satisfy the norm bound: 

sup lln^^^ll < r. 

k,(j 

Then B is well defined as a strong limit and is bounded in operator norm by 

\\B\\ < cr||?7||i. 

Proof. We check the conditions for Proposition IA.2I (applied for integrals rather than sums). 
Due to the intertwining relations between the unitaries rir^^ ^ and the multiplication operators 
gg^, we will then have a bound from above by an integral of multiplication operators. We must 
show that 1(^1 + G2) is bounded, where 

Gi= dk da |r| r^jqj:^J and G2 = dk da {q^^^^ r*jT^\. 
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The integrand of Gi is the multiphcation operator 
and the integrand of G2 is 

= \n^,^^{P + ak-k)M^[J+y^P,^',J + y^P)\ (3.2) 

where x^-^^ = Xj^o- + ycr(ao- — /), and we have used that TkM(P) = M{P — k)Tk- 

However since the operators in the integrand of Gi are all multiplication operators in P, 
bounding a sum on them in the operator norm can be computed as a supremum in the following 
way: 

IICill < sup II j dk j (i(T|ri^^(F)||r^(xi^^A; + y„F,X2,^A^+y^P)|||s(Hi) 
p J Jso„ 

< {sup \\n j: ^{P)\\ B{Hi)) { dk (ia |r7(xi,<^A; + y^P, xa.^A; + y<^F)|) (3.3) 

P ' P J JSO„ 

A similar result holds for 6*2 • Now applying Lemma [A. II to (13.31) along with our conditions 

on xi^o-, X2,o-, and n^^{P) we get the bound HGiH < rc||?7||i. 

□ 

Corollary 3.2. The integral of operators ( [^.5]) converges strongly to a hounded operator with 
norm less than or equal to IIpIIi- 



The bound in the above corollary in not sharp, since in Proposition (12. 2p we show that 
B = Tr2[pA^]. Thus ||P|| < ||p||i||Sa — /|| < 2||p||, since Sa is unitary. 

Proof. We apply Proposition (13.11) with n^^(P) = S'a(|A;|), rj = p, = a, xio- = 1 + A, 

X2_o- = / + cr, and y^- = A. |ng^(P)| < 1, so we can take r = 1. All determinants involved are of 
operators of the form ai + Xa2 where o"i, o"2 G SOn, so these determinants have a lower bound 
of (1 — A)". Hence we can take c = (1 — A)"*^. 

□ 

Proposition 3.3. Let G E BiHi ® L'^(W),nr L'^iW)), and ip : B{l-Li ® L'^{W),nr ® 
L2(M")) ^ B{n'^ L2(R"),7^o (g) L2(M")) has the form 

^{G) = J2[dk [ da,da2 Ul^^ h* ^ G g^^,^^^ f/,,,^ , 

where U^^ acts on the L^(M") tensor as Uj:^ = t^D\3^ t* ^, and h-^^ and gjj:^ are multiplica- 
tion operators in P of the form: 

hk,^ = ''fljP^'^f^^''^^-^ + ^2,<x^), <^nd g^,:^^ = nfl^{P)vf\^,J + X2,.P). 
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in BiTiul-L^) and B{Hr,'H^), respectively, and finally r]f\ rff^ G L^(M"). We will require that 



In the above, xi^o-, X2,o-, a^- G M„(R), bo- G G'Ln(M), the family of operators n . l and n . l lie 



inf I det(xi,o + ^2A^a^^<y - ^))l > - 



a 



and sup^.g^^ ll^-'ill'Sup,,fcV \\n^'^J < r. 

In this case, the integral of operator converges strongly to an operator ^{G) with the norm 
hound 

MG)\\<cr'\{\\T4, + \\T4{)\\G\\. 
where = \r]f'^){'>]j^\ for e = 1, 2 and || ■ ||i is the trace norm. 

Proof. We work towards showing the conditions of Proposition IA.3I with sums replaced an 
integral-sum. We thus need to find bounds for the operator norms of 



J2[ da I dkUlJgll^^^^U^,^^ ^nd J2 [ da! dkUlJh\l^^^^U^,^^. 

, JSO„ JK" . J. . JSOn "'K" 



(3.4) 



Is*!^^ = l^l^fc ^ multiplication operator with elements in B('H'i,'H'i) or an element in 

B{L'^(W^) ^Hi, L^(M") ^Tii). Conjugating with we get only multiplication operators back: 

With the calculations for bounding integrals of multiplication operators as in the proof of (13.11) . 
we get the bound 



{G)\\ < \cr'Y.^hT\\l + ll^f ll2)l|G|| = \cr\mV + r2||i)||G||. 



□ 



Corollary 3.4. The integral of operators ^2. (A) converges strongly to a limit with operator norm 
hounded hy ||p||i ||G|| (^)". 



Proof. We apply Proposition [331 in the case where = cr, bo- = xi^o- = f^, ^2,0 = '^^j+x^ 



In this case \n^'^l {P)\ and \n^'^l {P)\ < (1 — A) 2 so we can take r = 1. Also Xi o+xo a-(b„"'^ao — 

/) = ^i±^ and ||(^ij^)"i < 1, and hence | det(^it^)| > (1)" = 1 independent of A and a, 
so we can take c = 1. Hence by (13. 3p . we have our conclusion with a bound HpHiHGIKl — A)~". 

□ 
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4 Reduced Born approximation with third-order error 



In this section, we will prove Theorem 14. 2[ To make mathematical expression more compact it 
will be helpful to have the dictionary below. In the following expressions A, r G M"*", a G 5'0„, 
and A;, P G M". 



Dictionary of vectors in Dictionary of matrices in M„(IR) 

Q ,7_, ( p\ _ ^(l+A)-A-r(a-/) r X^«T-/p q „ _ (l+A)(r+(l-r)a) 

^k,a,r,Xy^) - iTA ^^T+A^ ^3,(7,r,A - <^(i+A)-Ar(<7-/) 

4- 4:,a(^) = TTa^-TTa^ 

Now we will list some relations between the vectors. The significance of these relations will 
become apparent once we begin doing calculations. 

Relations 

Rl- ^+^=TT^%,A(^) + TiA^ 

R3. k + P = Ci^a,r,X%^^^^^iP) + C2,a,r,xP 

In the proof of fl4.2p the analysis is organized around the fact that certain expressions are 
bounded. In the limit A — 0, expressions of the type ^5'a(-) will be a source of unboundedness, 
and p and G will have to be constrained in such a way as to compensate for this. The following 
expressions, defined for dimensions n = 1,3, are uniformly bounded in P,k G M, a G {+,—}, 
< r < 1, and < A: 

Ei{P,k,r,X) = ^-^^ jSx{\k\), (4.1) 

E,{P,k,a,r,X) = ^gA(Mg,,(P)|), (4.2) 

EsiPXX) = ^^"""^'^'g"'^"' y^A(|crg,,(P)|). (4.3) 

Their boundedness can be seen by using (Rl) to rewrite k in terms of a^^, ^(P) and P for 

Ei(P,k,r,X), (R2) to write d^^{P) in terms of Vk,a,r,x{P) and P for E2{P,k, a,r, X), and for 

E'i{P, k, A), dj: ^{P) exphcitly defined in terms of k and P. 

A second-order Taylor expansion of the scattering coefficients gives: 
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Dim-1 



Dim-2 



= A(l - X)i7rl - i\H\^ + ln(^)) - ^tt (4.4) 

Dim-3 

= ~ -A(l - A)2^//c - 2\H^k'' 

— — I ~r 

We can summarize the above expressions as 

Sx{k) ~ - A)c„fc"-2 - y c^A;2("-') - 5„,2^A2/'(7 + ln(|)), 

where Ci = a, C2 = tt/, and C3 = 2/. We will use the following simple lemma. 
Lemma 4.1. Let k,K eR and k,K e M^. 

1. We have the inequality 



< 2^ — < 



+ a 



^(A; - XK)^ + ^A2 



^. anc? /or dimension one the scattering coefficient satisfies 
5. anc? 

^. /or dimension three, the scattering coefficient satisfies 
k — \K ^ 4/ 

I'^'^Itt^I^ " (-2i/A|/c|)| < A^(i^(i + ^i^Dd^i + m- 

Proof. (1) follows by evaluating the critical points in A. (2) and (3) follow with an application 
of (1). 

□ 
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Define the following weighted trace norm || ■ \\^tn for the density matrices on the single 
reservoir particle Hilbert space p: 

||p|Utn=||p||l + $^ J2 \\\P^'^'[X^,[X„p]]h 
€ l<ij<n 

n 

+ ^^|||F|«-2+^X,pX,|Fp-2+^||,+ |||P|2(n-2)^|P|2(n-2)||^^ (4.5) 

£ j = l 

where the sums in e are over {0, 1} for dimension one and {—1, 0, 1} for dimension three. Notice 
the contrast between dimension n = 1 and n = 3 with respect to the weights applied in the 
norms for the absolute value of the momentum operators |P|. For n = 1, \\p\\wtn will blow up 
if p has non-zero density of momenta near momentum zero, while for n = 3, ||/9||«;tn can blow 
up if the momentum density does not decay fast enough for large momenta. This difference 
in requirements for different dimensions can be seen also in the formulas (11.5111.81) . The norm 
||/9||^t„ is not really asymmetric with respect to operators multiplying from the left and the 
right when p is self-adjoint. 

Theorem 4.2. Let e{G,X) be defined as in ( fi.ij) . then there exists a c s.t. for all density 
operators p E Bi{L\W)), G G B{L^{W)), and < X 

HG,X)\\<cX^\\pUn\\G\\^n. (4.6) 

Proof. We will prove the result for density operators p with a twice continuously differen- 
tiable integral kernel p{ki,k2) in the momentum representation, and a spectral decomposition 
P ~ of vectors fj{k) that are continuously differentiable in the momentum rep- 

resentation. Since such p are dense with respect to the || ■ \\wtn, the result extends to all p with 
llpll^to < oo. By f IB.ip . the Vi, V2, A operators and the map (p are well defined for all p with 
||/9||^f„ < 00 and they vary continuously as a function of p with respect to the norm || ■ \\u,tn- 

Our challenge is to expand the expressions we found in Propositions 12.21 and 12.31 in A, until 
we reach our second-order Taylor expansion while making sure that we only throw away terms 
which are bounded as in (14.61) . We will organize our analysis using the expressions (14. ip . (14.21) . 
and (14. 3p . in conjunction with Propositions 13.11 and 13.31 to effectively transfer the conditions 
for the boundedness of the differences in our expansions to conditions on p and G. Both of the 
expressions (12.50 and (12. 6p have multiple sources of A dependence. If we expand the expressions 
involving p and fj first for (12. 5p and (12.61) respectively, then the resulting expressions left to 
expand will be summable in the operator norm and thus not require the heavy preparation 
involved with the use of Propositions 13. Il and l3.3[ Breaking Tr2[/9S^G'SA] into parts and dividing 
by A we just need bound the differences 

jTi2[pAl]G - i^V^ + iXV2 + t^{A, P} - ^v^(/))G, and (4.7) 



jTT2[pAlGA]-X^{G), (4.8) 

where there is a similar expression to dUD for \GTT2[pAx\. We begin with (gTD, and will have 
to bound a sequence of intermediate differences. The main differences are the following: 
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Difference 1 



\\tt2[pAI]G- [ dk [ daT^,T;^^{pil,ak) + X{P + k)VTp{k,ak))\sxm)G\\, 



Difference 2 



[ dk I daTkT*^{p{k,ak) + X{P + k)VTp{k,ak)) 



(^5.(1^1) - ((1 - A)c„|^p-« + ^c^|^r-")))G||, 

Difference 3 

11/ dk daTf:T*^{p{k,ak) + X{P + k)VTp{k,ak)) 
Jr" Jso„ 

((1 - A)c„|fc|^-" + ^c^l^p^^-"))^ - {W^ + iW^ + ^{A F} - \p{^l))G\\. 
By the differentiability properties of tlie integral kernel p, 
p(fc + A(P + A?), aA? + A(P + A?)) = p(fc, ok) + A(P + k)\/Tpik, ak) 

+ X\P + kf' [ ds [ drVfp{k + X{P + k)r,ak + X{P + ky), 
Jo Jo 

where Vf g{x,y) is 2 tensor of derivatives with 

(\7^\(^ - hm ('^T9)^ix + hej,y + /tej) - (VrgUx, y) 

y\x,y)){ij) — ™i • 

The first difference can be rewritten as 

A^ fds fdr\\ I dk I rfar,^r;^^(l + |F|)(P+^)«^(5„3+|5^.^^_^(P)|n~2) 
Jo Jo Jr" J so„ 

Using (Rl) and expanding the tensor: (a^^^(P)+P)® a single term has the form ag^;^(P) P*^ ™. 
Note that the order of the tensors does not matter in this situation, since the whole vector is in 
an inner product with V*^ p, and partial derivatives commute. Now we apply Proposition 13.11 
with a single term: 



ri = {6n,3 + \kr')\kr{vfp),,,, 
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Finally with (O) we get the bound X'^C\\{6n,3 + |P|"-^)|Pr(V|V)i,fc||i|| I^P'^l^ + \P\)G\\, 
for some constant C. Note that VtP = ii^P — pX). 

The second difference can be bounded for dimension-one using the inequality 

A \k\ 2|fc|2" - |A;|3 

and for dimension three using the inequality 

I^S-Ad^l) - (2^(1 - A)|^| - 2A/|fcn| < 2\H^\k\''. 
At 

Finally, the last difference comes down to bounding the cross term: 

\\ f dk [ daT^T*^\iP + ak)VTp{k,ak){X^Cn\k\'''^ + X^cl\k\^^''-^^)G\\. 

Jr" JsOn " 



The bound for the above term follows from flA.ip and that / dk p{k ,k)^\^'^'^ = |||-P|" ^p|-P|" 



The ^G'Tr2[pAA] is similarly analyzed so now we study fl4.8l) . Again we have three main 
differences. There is a A dependence in tti^j:^^, Uj^^^, and S\{\dj: ^{P)\). It is most convenient 
to begin expanding mj^i:,a,x first. 



Difference 1 



|^Tr2[pA*GAA] -y, \ I dk f da, da2 U*^^ ^ det(/ + Aai)-^ f^{k) 



'SOnXSOn 



Difference 2 



II V f dk f da^da^Ul f,{a,k)\x-Uet{l + Xa,)-"^Sx{\d^,^{P)\)G 

■ JR" JsO„xSO„ ' 1' L 

S,{\^,{P)\) det(l + Xa^)--^ - cl \kf-''^ g] f,{a,k) U^.^^J 

Difference 3 

JSOnXSOn 



Using the differentiability of /j's 
f.{a,k-X'^^{k + P))= f^{a,k) + 



A (^) {k + P) drS/f, {a,k + rA (^) {k + P)) . 



(4.9) 
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The first difference 



||iTr2[pA^G'AA] -J^ [ [ dai da2 U^^ ^ det(l + Aai)-^ fj{k) 

is less than 

X^cl I dr|| V I dk I daida2Ui^ ^ det(l + Aai)"^ 

Jo j Jr" JsOnXSOr, 

- /DA , ^ fjiVk,cT-,,r,xiP)) 

(7— y)(c<xi,r,At^fc>,,,,A(^) + C2,<xi,r,AP))- 



«3 + 

m - (P) 

E^iP, k, au r, A)(l + |P|) G (1 + |F|)i?2(P, k, a^, r, A) J, ^ ^,^,,11 , 

where we have rearranged to substitute in the E2{P,k, a,r, X) expressions and used (R3) to 
rewrite k + P. Two apphcations of Proposition 13.31 corresponding to C(j-^^r,xVj:^^^{P) and 
C2,ai,r,xP will give us our bound. For the C(ji^r,\'Vj: ^ ^ ^{P) we use Proposition (13. 3p with 

vf\k) = {6n,s + \kr')\k\\VMk)\, 

Vf){k) 



= det(/+Aai)-5^2(P,A;,a,r,A)c,„,,A 



vf\k) = (5„,3 + |^r-')/.(fc), 



|V/,(A;)| 



nfl^^iP) = det(/+ Aa2)-5^2(P,fc,(T2,r,A), 
Hence the term is bounded by a constant muhiple of 

x'{j2mi<3Mpr')x,px,{6^^^^^^^ 

The C2^ai,r,xP term is bounded by a constant muhiple of 

x'{Y^\\{6r.,sMpr')x,px,{5^,,+\Pr^^^^^ 

The next intermediary difference has the form: 
III] / dk [ da,da2Ul^ det(l + Aai)-^/,(^)5,(|^%,(P)|) 

■ JR" JSOnXSOn ' ' 

Gifi^k,.,,xiP)) - fi^2k)) S,{\dk,,iP)\)det{l + Xa2)--- U^,^,J\. 
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Expanding .f{v^^^x{,P)) — f{o'2k) as in (14. 9p . we can apply a similar analysis to the above, 
except that for the left-hand side we organize around E^^P, k, A) rather than E^^P, k, a, r, A). 

Due to fj{aik)fj{a2k), the second difference is summable, and we do not need to prepare 
any more applications of Proposition 13.31 We begin by bounding 

||5^A/c^ [ dk [ da.da^UI det{l + Xa,)--^ 



1^1'"" cnJ .mn r^]fA<^2k) 



We observe the inequality 



Sx{\daP)\)-G 



\k\ 



det(l + AcT2)"5 



2-n ^ ^'"^a'^' 



|jL|2-n It-I2-" 
1 uAk 



S,{\d^^,{P)\)G^S,{\d^,,{P)\) -G]\\< 



2-n 



Cn. A Cn 



-5,(|c%_,(P)|)-^)G(|P| + /)E3(P,^,A)| 



1 2-n 



+ I|g(^^a(|c%,,(p)|) + ^)||. 

By (3) and (4) of (14.11) . the right-hand side is bounded by a sum of terms proportional to 
A|fc|''("-2)|||P|'iG'(/ + for r = 0,1,2, ei,e2 = 0,1. Bounding the above integral is then 

routine and requires that || |Pp'-'^~^^p|Pp*^"~^)||i. The last thing to do for the second difference 
is expanding det(l + A(Ti)~2 and | det(l + \a2)\^^, which does not pose much difficulty. 
For the third difference, we will need to work with the D i+\a term. 

l+A 

\\UUxGU^,^^^^,-T:^^^r^,Grlr^J\ < + \\GiD._^ - 1)1 

l + A -L + A 

since Uj:^^ ^, t^j^, and are unitary and D*+;^^^r^ = T i+x,y ■ Di+^ satisfies the integral 

' ' l+A ^+-^ l+A 1+-^ 

relation ^ 

Di+x.=I+[ ds {4- log( ^ ^ ^^ )Di+s.P^ Xj, 
Jo ds 1 + s ~ ^ 

and hence 

\\\{D._^ - I)G\\ < ( sup ;f log(^^))|| WP^^M- 
A 1+^ o<s<A as 1 + s ^ 

— *j 

The third difference is then bounded by a fixed constant multiple of A^|| iPl^^^plPl^^^Hi || 

□ 

APPENDIX 



A Hilbert spaces and operator inequalities 

Lemma A.l. Let rj be a trace class operator on MJ^ with continuous integral kernel r]{ki,k2), 
A, A' e GL„(M), and a, a' e M", then 

£ rff |r,(Af + a, A'f + a-) I < + 
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Proof. Let p = J2j then we have 

/ dx \p{Ax + a, A'x + a')\ = dx\^ Xj fj [Ax + a) gj [A'x + a') \ 
Jr" Jr" j 

where the inequahty follows from 2ab < + b'^ for a,b & M.'^ and completing the integration. 
The equality on the left-hand side of ( lA.ll) is formal for a general integral kernel p{ki,k2) 
which is defined only a.e. with respect to joint integration over ki, k2, but with our continuity 
condition it is well defined. 

□ 

Proposition A. 2. For n eN, let An G B{T-l) for a Hilbert space Ti. and 

n 

be weakly convergent to a bounded operator with norm c. Then An is strongly convergent 
to a bounded operator X with \\X\\ < c. 

Proof. Let g E H and with the polar decomposition [TT] An = Un\An\, then taking a tail sum 



\\^Ang\\2= sup {h\^Un\An\^An\'^g) 
JV il^l|2 = l ^ 

OO oo oo 

< sup 5^|||An|^f/n/^||2|||A„|^^||2< SUp ( (/iP: | /l) ) ^ ( ^ P„ | ^) ) ^ 

Il^ll2=l N \M2=1 n=l N 

where the first and second inqualities follow by two different applications of the Cauchy- 
Schwartz inequality. The right-hand side then tends to zero for large by our assumptions on 
the series ^ |^n| and l^nl- The operator norm bound can be seen from the same calculation 
with a sum over all n rather than a tail. 

□ 

Proposition A. 3. Let An, Bn for n E N be elements in B(H) for a Hilbert space H such that 
J2n ^n^n and Ylin ^riBn convcrgc weakly to bounded operators with norms less than c, then the 
sum 

^{G) = Y,KGBn 

n 

is strongly convergent to an operator with norm less than or equal to c\\G\\. 
The proof follows a similar pattern to that of Proposition IA.2I 



B Norm bounds for Vi^ V2, A and (p 

The following lemma shows that the limiting expressions f 1 1.3 1) vary continuously with respect 
to the density operator p in the || ■ topology. It allows the limiting expression to be defined 
for all p with ||p||«,tn < oo without the additional assumption that the integral kernel of p in 
the momentum representation in continuously differentiable. A somewhat weaker norm than 
II ■ Wwtn would suffice. 
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Lemma B.l. Let V\, V2, A, and ip be defined as in ( f j.5]) -( l77^) for ape i3i(]R"), p > 0, 
with continuously differentiate integral kernel in the momentum representation, then there is 
a constant c > such that for all p and j 

llFlll, 11^211, ll^ll, ||[P,,A,]||, lly^ll <c||p|Ui„. 

Proof. By an argument similar to flA.ip 

ii^iii < c^\\\p\p\\i, \\A\\ < cnY,\\\pr'Xjp\\u M < c^ii|p|"-vi^r~ii, 
m < cnj2m\pr\[Xj,p]}\\i, m,AM<cnm\pr',[x,,p]]\u, 

where we have used that \\p\\ = \\p{I)\\ since p is a. positive map. By p being self-adjoint, we 
have inequalities such as 

ii^.i^r"V^.iii < ^(ii^.p^.iii + \\p,\pr'p\pr'p,h), 

and then that Pj < \P\. Finally, since p is positive dkp{k,k)\k\'^^ = || iFl^plFl"*!!!, and 
inequalities of the form 

iiip|>|priii<iipiii + iiipivi^nii 

follow, where r,s have the same sign and \r\ < \s\. Hence ||p||iotn bounds the expressions for 
Vi,V2, A, [Pj,Aj] and p. 

□ 
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